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Abstract 

First,  a consistent  perspective  for  the  formulation  of  constitutive  dyadics  for  biaxial  medi- 
ums — for  the  anisotropic  dielectric  and  the  full  bianisotropic  cases  — is  provided.  Then, 
the  connection  between  the  existence  of  closed-form,  infinite-medium,  dyadic  Green  func- 
tions and  the  factorization  properties  of  certain  scalar  differential  operators  is  explored  by 
focusing  on  a special  type  of  homogeneous,  anisotropic,  dielectric  medium.  Its  anisotropy  is 
of  a higher  degree  of  complexity  than  an  uniaxial  medium’s  but  falls  short  of  a fully  biaxial 
medium’s. 

1.  From  Uniaxiality  to  Biaxiality 

An  anisotropic  dielectric  medium  of  the  simplest  type  has  a relative  permittivity  dyadic  e that 
is  uniaxial.  Stated  as1 

guni  = eal+ebuu , (1) 

it  employs  two  complex- valued  scalars  (i.e.,  ea  and  eb)  and  one  unit  vector  (i.e.,  u)  which  is 
parallel  to  the  sole  distinguished  axis  of  the  medium  [1].  Extending  the  structure  of  (1)  to  the 
permeability  dyadic  and  the  magnetoelectric  dyadics  as  well,  we  arrive  at  a uniaxial  bianisotropic 
medium  [2]. 

Generalization  from  medium  uniaxiality  to  biaxiality  requires  the  introduction  of  a second 
distinguished  axis.  Recently,  we  [3]  put  forward  a consistent  approach  to  that  issue  for  biaxial 
bianisotropic  mediums , delineating  their  frequency-dependent  constitutive  relations  as 


D(x)  = e0eb.  • E(x)  + (^e0/p0)abi  * B(x) , (2) 

H(x)  = (\/eo/Mo) * E(x)  + (1  /po)  • B(x) . (3) 

Here,  the  four  constitutive  dyadics  are  given  by 

gbi  = Ca  I + eft  («m  un  + Un  Um)  , (4) 

abi  =aal  + ab  (um  un  + un  um)  , (5) 


Hn  this  paper,  vectors  axe  in  boldface;  dyadics  axe  underlined  twice;  a • b = aibi ; A = ab  is  a dyad 

with  elements  Aij  = ca  bj\  I is  the  unit  dyadic  and  0 is  the  null  dyadic;  the  superscript  indicates  inversion  of 
dyadics  and  differential  operators;  while  x and  x'  are  the  observation  and  the  source  points,  respectively. 
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P..  = Pal  + Pb  (UmU„  + UnUm)  , (6) 

= 01  — 

X..  = XaI  + Xb  (umUn  + U„Um)  , (7) 

= 01  “ 

with  um  and  un  as  two  unit  vectors  that  are,  in  general,  neither  parallel  nor  anti-parallel  to 
each  other.  The  simplification  um  = ±un  leads  us  back  to  uniaxial  bianisotropic  mediums. 

The  two  unit  vectors  in  (4)-(7),  and  both  corresponding  distinguished  axes  [1],  are  common 
to  all  four  constitutive  dyadics.  Furthermore,  these  dyadics  contain  eight  complex-valued  par 
rameters:  ea,  eb,  aa,  ab,  pa,  Pb,  Xa  and  Xb,  while  the  angle  <p  = cos-1(um  • un)  is  real-valued. 
Therefore,  any  biaxial  bianisotropic  medium  is  characterized  by  8 x 2 + 1 = 17  real-valued 
constitutive  scalars.  A redundancy  in  this  scheme  is  filtered  out  by  the  Post  constraint  [4]: 

Trace  (o_bi  - P^j  =0  =>  3 (a0  - Pa)  + 2 (ab  - pb ) cos  ^ = 0 . (8) 

Hence,  (4)-(7)  actually  contain  just  15  real-valued  independent  constitutive  scalars.  Further- 
more, as  presently  defined,  the  biaxial  bianisotropic  medium  is  a nonreciprocal  medium. 

A biaxial  anisotropic  medium  [1]  arises  as  a special  case  when  = P_  = 0 in  (4)-(7).  The 
biaxiality  is  thus  purely  in  the  dielectric-magnetic  properties,  and  is  described  by  4 complex- 
valued constitutive  parameters  plus  the  real-valued  angle  <p.  As  the  Post  constraint  (8)  is  then 
trivially  fulfilled,  a biaxial  anisotropic  medium  is  uniquely  described  by  9 real-valued  constitutive 
scalars. 

Perhaps  the  most  attractive  feature  of  the  representation  (4)-(7)  is  found  in  the  possibility 
that  one  and  the  same  orthogonal  transformation  is  able  to  diagonalize  all  four  constitutive 
dyadics.  Full  details  about  that  feature,  as  well  as  a comprehensive  study  of  electromagnetic 
wave  propagation  in  biaxial  bianisotropic  mediums,  are  available  elsewhere  [3]. 

2.  Dyadic  Green  Functions 
The  dyadic  differential  equation 


£(V)  • £(x,  x')  = <5(x  - x')  I , (9) 

constitutes  the  standard  field  problem  for  the  electric  field  phasor,  with  the  exp(— iu)t)  time- 
dependence  implicit  throughout.  The  dyadic  Green  function  (specifically,  of  the  electric  type)  is 
denoted  by  £(x,  x'),  while  S(x  — x')  is  the  Dirac  delta  function.  The  specific  form  of  the  dyadic, 
second-order,  differential  operator  L(V)  depends  on  the  type  of  medium  being  considered;  see 
Ref.  [5]  for  the  relevant  exposition  of  the  general  Green  function  technique. 

For  brevity’s  sake,  we  limit  our  attention  here  to  anisotropic  dielectric  mediums,  whose  per- 
meability dyadic  equals  X and  whose  magnetoelectric  dyadics  equal  Q.  A closed-form  expression2 
for  Q(x,x')  for  an  uniaxial  dielectric  medium  is  available  in  textbooks  [1],  but  none  is  known 
to  exist  for  a biaxial  dielectric  medium.  The  latter  is  not  surprising  in  view  of  some  results 
pertaining  to  the  so-called  determinant  operator  of  L(V). 

An  important  step  towards  a closed-form  solution  is  to  determine  whether  the  determinant 
operator 

B*,  = i(V)  • £^(V)  = ia4j(V)  • £(V) , (10) 

can  be  factorized  into  a product  of  two  second-order  operators,  the  adjoint  operator  £a(^(V)  in 
(10)  being  always  of  the  fourth  order  by  virtue  of  the  structure  of  the  Maxwell  equations.  It  was 

2 The  dyadic  Green  function  should  be  expressible  through  simple  mathematical  functions  which  will  most  often 
be  scalar  Green  functions  of  second-order  Helmholtzian  operators  and  derivatives  as  well  as  linear  combinations 
thereof.  It  does  not  include  representations  in  terms  of  integrals,  because  such  representations  can  be  trivially 
achieved  with  spatial  Fourier  transformations  due  to  the  linearity  of  (9). 


5 


first  stated  in  Ref.  [6]  — see  also  Ref.  [7]  — that  the  determinant  operator  can  be  factorized 
only  if  the  relative  permittivity  dyadic  has  the  structure 

s,„<j  = Ai+ab-  (“) 

where  A is  a scalar  while  a and  b are  vectors.  Parenthetically,  the  formulas  given  in  Refs.  [6,  7] 
pertain  to  the  more  general  anisotropic  dielectric-magnetic  mediums. 

The  form  of  e^act  stated  above  is  only  sufficient  but  not  necessary  for  factorization.  Fur- 
thermore, the  relation  between  factorization  and  the  availability  of  closed-form  solutions  is  not 
clear.  In  fact,  closed-form  expressions  for  dyadic  Green  functions  could  not  be  obtained,  despite 
factorization,  for  certain  types  of  uniaxial  bianisotropic  mediums,  as  was  first  observed  in  Ref. 
[8]  and  explored  further  in  Refs.  [9,  2]. 

3.  More  than  Uniaxial  — Not  quite  Biaxial 

We  now  consider  an  anisotropic  dielectric  medium  whose  relative  permittivity  dyadic  is 

g = CoI+eftUmUn,  (12) 

where  um  and  un  are  distinct  unit  vectors.  The  right  side  of  (12)  is  clearly  equivalent  in  form  to 
£fact-  In  general,  the  chosen  medium  is  not  reciprocal  as  e does  not  equal  its  transposed  dyadic. 
The  uniaxial  medium  defined  through  (1)  appears  as  a specialization  of  (12)  when  um  = un. 
However,  (12)  also  has  a connection  to  a biaxial  structure:  Upon  decomposition  of  its  right  side 
into  symmetric  and  skew-symmetric  parts,  (12)  can  be  rewritten  as 

| = ea  l + ^ (um  un  + un  um)  + j (um  un  - u„  um)  . (13) 

The  first  two  terms  on  the  right  side  of  (13)  have  the  exact  biaxial  structure  of  (4),  whereas 
the  last  term  has  a typical  gyrotropic  form.  Nevertheless,  neither  the  biaxial  dielectric  nor  the 
gyrotropic  dielectric  medium  can  be  obtained  from  (13)  as  special  cases  — because  the  last  two 
terms  on  the  right  side  of  (13)  are  far  too  intimately  linked  to  allow  those  specializations. 

For  the  medium  characterized  by  (12),  we  obtain  [10] 

Ladj(V)  = #m£e(V)  - k2r  (V  x un)  (V  x um)  , (14) 

Hdet  = -k2  (1  + r um  • un)  He  Hm  , (15) 

where  k 2 = uj2e0eap0)  the  ratio  r = eft/e0,  and  the  dyadic  operator 

Le(V)  = V V + k2  (1  + r um  • u„)  (l-  — — um un)  . (16) 

V 1 “h  T Um  • Un  / 

Of  the  two  scalar,  second-order  operators  appearing  in  (14)  and  (15),  Hm  — V2  + k2  is  a 
standard  Helmholtz  operator  due  to  the  magnetic  isotropy  of  the  medium,  but 

tfe  = V2-— — (V  x um)  • (V  x un)  + k2 , (17) 

1 + T Um  * Un 

is  only  a Helmholtzian  operator  as  it  reflects  the  dielectric  anisotropy  of  the  medium.  As 
anticipated  for  the  chosen  medium,  it  is  apparent  from  (15)  that  Hdet  indeed  factorizes  as  a 
product  of  two  second-order  operators  [10]. 

Farther  manipulations  then  lead  to  the  complete  representation  of  £(x,  x')  in  the  form 

kHl+r^-u „)  [i,(V)9«(x,x')  + iVM(x,x')]  , 


£(x,x')  = 


(18) 
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where  the  scalar  Green  function 

5e(x,x') 


1 exp  [ikD(x,  x')] 
y/axayaz  47rL>(x,x') 


(19) 


satisfies  the  differential  equation  Hege(x,xr)  — -<5(x  - x'),  and  the  structure  of  the  modified 
distance  function  D(x,x')  (which  also  contains  aXiy}Z)  is  detailed  in  Ref.  [10]. 

What  remains  unknown  in  (18)  is  the  dyadic  function  for  which  the  differential 

equation 

HeHm  M(x,  x')  = (V  x un)  (V  x um)  S(x  - x') , (20) 

must  be  solved.  As  discussed  elsewhere  [10],  no  closed-form  solution  for  M(x,  x')  appears  to 
emerge  from  (20). 

The  medium  considered  here  is  the  most  general  type  of  anisotropic  dielectric  medium  for 
which  the  determinant  operator  is  the  product  of  two  scalar,  second-order  differential  operators 
of  the  Helmholtzian  type;  yet  no  closed-form  expression  for  £(x,  x')  appears  to  exist  for  this 
medium!  Thus,  we  conclude  that 

(i)  while  all  currently  known  closed-form,  infinite-medium  dyadic  Green  functions  are  based 
on  factorizable  determinant  operators,  factorization  is  not  a sufficient  condition  for  the  existence 
of  closed-form,  infinite-medium,  dyadic  Green  functions;  and 

(ii)  within  the  class  of  anisotropic  dielectric  mediums,  the  uniaxial  dielectric  medium  (or  any 
medium  that  can  be  reduced  to  such  a medium  by,  for  example,  affine  transformations)  remains 
the  most  general  medium  for  which  a closed-form,  infinite-medium,  dyadic  Green  function  has 
been  derived.3 
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